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1 Introduction 



Matrix model provides a new paradigm for thinking about fundamental theories of physics. It 
originates in the observation H that the massless modes propagating along the world volume 
of N coincident D-branes are those of the supersymmetric Yang-Mills theory, obtained by the 
dimensional reductions of the D = 10 N = 1 theory down to p + 1 spacetime dimensions. 

According to so-called Matrix conjecture [||, + 1-dimensional reduction can be regarded 
as the discrete light cone quantization of M-theory in which the spacetime is compactified on 
an almost light like circle. This proposes a concrete, a nonperturbative definition of quantum 
gravity, and quite remarkably, the conjecture has found quite nontrivial support j|, [I], || |6| . 

Meanwhile, type IIB matrix model proposed by |7| plays somewhat complementary role. Its 
action is + 0-dimensional reduction of large TV super Yang-Mills theory in ten dimensions. The 
authors of || proposed a very interesting program to study dynamical formation of space-time 
using the type IIB matrix model. 

Despite much remarkable success of Matrix model approach, the question "Why and how 
such a simple model could describe our real world?" is still elusive. The main difficulty consists 
in the absence of built-in rules concerning "How to take large N limit." For example, in the 
case of Matrix model approach to 2d gravity there is a critical point g c , and continuum limit 
is possible keeping certain relation between N and g — g c (double scaling limit). For type IIB 
matrix model, however, the coupling constant g can be absorbed into the rescaling of the fields 
(at least classically) and there is no nontrivial fixed point. 

The result of matrix integration is just a number as it stands. To extract physical intuition, 
we need to separate field variables into two types: the classical background and the quantum 
fluctuation. 

In the spirit of Born-Oppenheimer approximation, the effective dynamics of the slow vari- 
ables (classical background) are of primary concern which is obtained only after fast variables 
(quantum fluctuations) are integrated out. This is the approach taken by many works. 

In Matrix models, however, somewhat different approach might be of considerable interest. 
Recall that in the usual analysis of quantum field theory, gravitational effects are almost always 
ignored, although gravitons are massless and never decouple. Gravitational degree of freedom 
are not integrated over, but regarded as fixed, classical background. This treatment is justified 
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simply because the dimensionful coupling is so small in the energy scale accessible by the current 
technology. Similarly for the observers living on the branes, natural time scale is set by that 
of quantum fluctuations rather than the dynamical time scale of the background = spacetime. 
Put it differently, "motion of the background is too slow to be treated quantum mechanically." 

In this paper, we will study the quantum dynamics of the Matrix model from the latter point 
of view, hoping our work provide some insight about how to take large N limit. The paper is 
organized as follows. 

In section two, starting from + O-dimensional matrix action, we derive fatgraph Feynman 
rules for the quantum fluctuations treating general multi D-instanton configuration as a fixed 
background. The usage of the Feynman rules is shown with an example. Although we will work 
in D-instanton backgrounds of IIB matrix model, we expect our analysis shed some light on 
general D-p branes in other Matrix theories as well, since type IIB matrix model compactified 
on S 1 is equivalent to the 1 + O-dimensional Matrix model Ji~0| ]. 

The matrix Feynman rules are very close to those in the usual ci-dimensional SYM. In 
section three, we will study a special backgrounds where D-instantons are concentrated along 
(i-dimensional sheet in the original D-dimensional spacetime. We will see finite N theories can 
be thought of as UV regulated versions of flat space Yang-Mills theory in which removing the 
cutoff is equivalent to letting N go to infinity. The crucial observation of this paper is that 
from Yang-Mills perturbation point of view, going to Matrix model can be thought of as a 
discretization of a momentum space rather than a coordinate space. This is shown explicitly by 
comparing Feynman rules. This is yet another manifestation of spacetime uncertainty [0] or 
UV/IR correspondence JT^]. 

For the ci-dimensional quantum field theory embedded in the Matrix model, the only source 
of divergence is the large N limit. In section four, we will study the structure of large N 
divergences in Matrix theory and relate it to the renormalizability of QFT in the usual sense 
of the term. We hope this line of argument give us a hint to deduce realistic physics from the 
Matrix models. Finally, we conclude with a discussion of our results and some implications. 
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2 Matrix perturbation theory around D-instanton background 

In this section, we elaborate the perturbation theory of the Matrix model around D-instanton 
background. Explicit forms of Feynman rules are derived and boson self energy diagrams are 
computed at one loop as an example. 

2.1 The type IIB matrix model and its gauge fixing 

Our starting point is the Euclidean type IIB matrix model, whose action is given by 

where X^ and ip are D-dimensional vector and Majorana-Weyl spinor respectively, taking values 
in N x N hermitian matrices.^ Throughout this paper, "Tr" denotes the trace taken over N x N 
matrix indices. 

The action enjoys the following symmetries 

• rotation invariance 

SX* = uj^X u , 5iP = iw'T^, {u^ = -uj Ufl ) 

• translation invariance 

SX* = (2.2) 

• M = 2 SUSY 

5^^= l -[X^X u ]T^e 1 , S^X lt = ie 1 T^, <*< 2 ty = e 2 , 5^X^ = 0. 

• U(N) Gauge invariance 

X» i ^ U^X^U, %b i ^ U-^ijPU, (U € U(N)) (2.3) 

• scaling property : 

X^XX", g^X 2 g (AeR >0 ) (2.4) 

1 We choose D = 10 type IIB matrix model just for definiteness. We could start from any model reduced from 
D-dimensional SYM. 
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As stated in Introduction, we decompose X^ as a sum of classical background part and 
quantum fluctuation part X^. X^ will be treated as fixed, classical number and we will be 
interested in the quantum field theory in this background, (ip is assumed to have no classical 
vacuum expectation value.) 

The background X^ must be a solution to the equation of motion, [X^, [X U ,X^]] = 0. We 
will consider the cases where all the X^'s are simultaneously diagonalizable by the gauge action 



( 



x» = x^ + x^, 



x» 



\ 



(2.5) 



X N ) 



The combination of D eigenvalues Xi = (xj, ■ ■ ■ xf) £ ~M, D is interpreted as the location of the i- 
th D-instanton. For a generic background where all D-instantons are separated from each other, 
all the symmetries listed above are explicitly broken. In particular, U(N) gauge symmetry is 
broken down to U(1) N and half of the M = 2 SUSY (generated by 6^') survives indicating the 
BPS nature of the background fl2,5|). 

Let us make a brief comment on the charges of the fields. All the quantum fields X^tpij, Cij, 
bij have a charge 



3^ 

0-10 



(2.6) 



( ■•■ 1 ■•■ -1 ••• 0) 

with respect to the unbroken U(1) N gauge symmetry. These fluctuations correspond to the open 
string stretching between D-instantons i and j. In particular, diagonal components X^^ipn^cu 
and ha are neutral. In fact, as we will soon see, their kinetic terms vanishe indicating they should 
be treated as collective coordinates rather than quantum variables, and thus need a separate 
treatment. Since these diagonal components could be absorbed into the shift of the background 
D-instanton configuration, incorporating these fluctuations would inevitably lead to the integral 
over the collective coordinates, which is beyond the scope of this paper. 

We will study the quantum theory of fluctuations as parameterized by the classical back- 
ground i.e. D-instanton positions {xi]f =1 . Plugging ( |2.5[ ) into the action ( |2.1| ) and using the 
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relation [X»,X»] ik = (xf - we have 

+ Yl^ Xik + X 3kY X ij XU jk X k-i\ + 7^2 E i X ij X jk X kl X li ~ X ij X jk X kl X li} 

i,j,k 

~h EWX-^} + i E^ rA ^**« - ^r^^}. (2.7) 

It should be noted that Xj's always appear as difference x% — Xj due to translational invariance 
( |2.2| ). Hereafter the notation x^- = x^ — Xj will be used to simplify the formulas. 

To setup a perturbation theory, convenient to work with the background field gauge. This 
is achieved by adding the gauge fixing term 

S 9.f- = -^Tr[lMf = —D^AJi} (2.8) 

accompanied with the Faddeev-Popov ghost term 

Sf.p. = ~Tr[X»,b][X' l ,c] 

g 2 

= 3 £{(*&) a *y<*} + 3 " (2-9) 

i,j g i,j,k 



The gauge fixing term fl2.Sj ) implies we have chosen a gauge such that the fluctuation Xjj is 
transverse to the relative vector Xij, i.e. Ylu x ij X ij = 
The gauge fixed total action is given by 

Stotai = S + S g j, + Sf.p. (2-10) 



Note that S g _f. in ( |2 . 1Q| ) cancels with the second term of (2.7). The perturbation is valid when 
D-instanton separation Xy is much larger than g 1 ^ 2 . 

2.2 Feynman rules 



Now we will derive Feynman rules from the gauge fixed action ( 2.1 Q| ) . Hereafter we will rescale 



the quantum fluctuations X^ — > gX-j, etc, in order that the coupling g is removed from the 
propagators and moved to the interaction vertices. 

The perturbative structure of large N gauge theories are naturally described in terms of 



double line representation [13| of Feynman diagrams, so called fatgraphs. One considers, as in 



i — i ^ 

fermion a ~+ (3 ^i'd'H~ 'J—^i. 



ij 

i 
3 



?host c-ijbji^ — 

............. x ^ 



Figure 1: Propagators 

Figs. |l] and |2[ a graph with the lines thickened slightly into bands which meet smoothly at the 
vertices so as to form an oriented Riemann surface with boundary. In this representation, the 
fields in the adjoint representation of U(N) are denoted by double lines. 

Each edge carries a label € {1, . . . ,N} corresponding to the basis of fundamental 

representation of U(N) (or its conjugate depending on its orientation). Here it is nothing but 
the label of a D-instanton. From the quadratic part of the total action ( p. 10 ), we can easily read 



off the fatgraph propagators as depicted in Fig. [T|. Note that the denominator is the squared 
distance between the two D-instantons connected by the fields. 

As for vertices, there is a crucial difference between ordinary (particle theory) vertex factor 
and fatgraph counterpart. In the former, interaction vertex of order k is invariant under all 
possible k\ permutations of lines, while in the latter, it is invariant only under k cyclic permuta- 
tions. Thus, for example, two Yukawa coupling diagrams in Fig. ^ should be distinguished from 
each other. 

In deriving the vertex factors, terms must be organized so that the index contraction should 
have manifest cyclic invariance. For instance, in order to deduce three-point and four-point 
vertex for X's, we must rewrite corresponding terms in (2/7) as follows 



^ \ " \ " y-Ml W»2 yfJ-S 
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2g 2 

i,j,k,i 



' \ / yV yV yv yv yfi yv yH yV\ 
^2 Z-> l A ij A jfc A W A & - A i, A jfc A j« A fi/ 

' y-Mi y^A 

72 2^ 2^ A £i A ij A jfc A M 
i,j,k,£ pi,— ,H4 



Overall factors 1/3 and 1/4 are cancelled by the cyclic symmetry of the vertices. 
Similar computation leads to the vertex factors listed in Fig. |2[ 

2.3 Example: one loop boson self energy 

In order to illustrate how matrix perturbation theory works, let us calculate the one loop contri- 
bution to two point function^ (X^X^) using the Feynman rules just derived. Relevant fatgraphs 
are shown in Fig. ||. In addition to the "external" D-instantons % and j, we need to incorporate 
an "internal" D-instanton k as in Fig. f|. 

Let us begin with the diagram (a). It can be seen as representing the history of an open 
string; the string ij splits into two pieces ik and kj, and reconnect in the end. Applying the 
Feynman rules, this diagram contributes 

a 2 1 

V k k\ X ik X jk 

x{6^( Xij - x 3k ) x + 5^{x ik - x kj y + 5^( Xjl - x ik ) K } 

x{5™(x kj - Xji ) x + S vX ( Xji - x tk ) K + 5 XK (x ik - x kj f}. (2.11) 

Here we denote by A and k the SO(D) vector indices associated to the upper and lower internal 
propagators respectively. Similarly, from the diagrams (b) and (c), we have 

(b) + (c) = -4 £( 4- + 4-) (1 - D)8, v . (2.12) 

x ij k yx ik x jk y 



Ghosts contribute 



U k X ik X jk 



2 no sum is taken for indices i or j; it is invariant under U(l) gauge symmetry. 
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i 




Figure 4: D-instanton configuration associated with the one-loop processes in Fig. ^ 
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Finally fermion loop diagram gives 

(D + (g) = £ ^ + ^;^"^ (2.i4) 

lX %3 k X jk X ik 

where dr is the size of gamma matrices in Z)-dimensions. 

3 Correspondence to quantum field theories 



From the sample calculation given in section 2.3, we notice a strong similarity between the 



matrix perturbation theory and a usual <i-dimensional QFT. Roughly speaking, relative brane 
position Xij corresponds to a momentum p whereas the sum over branes Ylk l°°ks like a loop 
integral f d d p. 

In this section, we will make this analogy more precise. In particular, we illustrate how d- 
dimensional gauge theories can be recovered from Matrix model, when D-instanton configuration 
has ci-dimensional flat directions. In this context, the flat directions should be thought of as 
momentum coordinates rather than spatial coordinates, contrary to naive expectations. This 
can be considered as an example of IR/UV correspondence [Q|. 

Infinite sums or integrals will pose a delicate problem of large N divergences. We will 
postpone discussing this issues to section |||. 

3.1 D-instanton distribution as discretized momentum space 

We begin with recalling a general structure of ci-dimensional gauge theory amplitudes. Let A^(x) 
be U(n) gauge fields represented asnxn matrices and A fM (p) their Fourier transform. Hereafter 
the symbol "tr" will denote the trace over nxn matrix indices. Any gauge invariant correlation 
function can be decomposed as a sum of basic correlation functions of the form 

(tr A^(jp l )A^{p 2 )---A^(jp k )) 1 (p 1 + ... +Pk = 0) (3.1) 

which is invariant under the cyclic permutation of momenta and Lorenz indices: 

P\^V2^ >Pk ->■ Pi, 

(3.2) 

Mi -* M2 -> > Ml- 

In the amplitude ( |3.1| ), the cyclic order of {(pi,Mi)}i=i has a definite meaning because U(n) 
indices are implicitly contracted. It is easy to check that there is one-to-one correspondence 
among the following three data: 
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P3,V3 



Pk,Vk 





Pi, Pi 



(a) 



(b) 



'■2 



Figure 5: Correspondence of a correlation function between gauge theory and Matrix model, 
(a): closed path formed by external momenta, (b): corresponding D-instanton configuration. 



(A) gauge invariant amplitude ( |3.lD , 

(B) ordered set of momenta and Lorentz indices {{pi, Mi), (p2, ^2), ■ ■ ■ , {Pk-> Mfc)} modulo cyclic 
permutation (3.2), 



(C) closed oriented path in M. d with edge labeling Ml, * * ' , Mfc (Fig- H (a))- 

We now come back to Matrix model correlation functions. In generic background ( |2.5[) we 
still have unbroken U(1) N gauge symmetry. Therefore, only gauge invariant "Wilson loops" 
such as 



\ ll«2 *2«3 *fcH' 



(3-3) 



can be nonzero. Actually, we can draw a corresponding loop as in Fig. || (b): two D-instantons 
located at Xi r and Xi r+1 are connected by a field (or an open string) j . The amplitude 

is invariant under the cyclic permutation of D-instanton positions and Lorenz indices: 



Ml M2 



-» Mfc -» Mi- 



(3.4) 



The correspondence between the two figures Fig. || (a) and (b) is now obvious; the momentum 
p r in gauge theory is identified with the D-instanton separation Xi r+l — xi r . In this way, we can 
add two new entries to the previous list of one-to-one correspondence: 



(D) D-instanton correlation functions (|3.3| ) 
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(E) loop passing through k D-instantons in M. d with labeled edges Hi,. ■ ■ ,Hk 
(Fig. | (b)). 

Clearly, arbitrary sequence of (i-dimensional momenta {pi, ... ,pk) can be realized as a loop in D- 
instanton configuration space, provided D-instantons are densely distributed over d-dimensional 
Euclidean space M. d . Conversely, if we start from finite N matrix theory, only discrete momentum 
points are available on gauge theory side. Of course, the error becomes smaller as the number 
of D-instantons is increased. Guided by these observation, we propose to identify a D-instanton 
distribution of Matrix theory as a discretization of momentum space seen by a Yang-Mills theory. 

Precisely speaking, the "momentum path" (Fig. || (a)) can determine the "D-instanton path" 
(Fig. H (b)) only up to overall translation, Xj — > Xi + c. This ambiguity can be resolved by fixing, 
say, their center of mass at the origin. 

So far we neglected the problem how non-Abelian U(n) gauge symmetry can be recovered 
from D-instanton picture. This will be discussed in section |3,4| . 

In sum, we have argued that when the background D-instantons are continuously distributed 
along M. d , there is a one-to-one correspondence between the gauge theory amplitude (pO) and 



3.2 Correspondence of Feynman diagrams 

In perturbation theory, both amplitudes Q3.1|) and (3.3) are expressed as a sum over fat graphs 



with fixed external lines. We now want to show that two computations, one as a Yang-Mills 
theory and the other as a Matrix theory, actually coincide for every fatgraph. To do this, we 
need to check the correspondence at the level of propagators and vertices. 

Consider a fatgraph T made of several propagators and vertices. Recall that the graph form 
an oriented Riemann surface with boundaries. Thus for a given propagator with an orientation, 
it is meaningful to talk about its "left-" and "right-" edges. 

Pick up a propagator and let i and j be its labels on left- and right-edges, repectively. In 
Matrix picture, the propgator represents a fluctuation Xf- connecting two D-instantons i and j. 



As in section 3.1, we identify the relative separation of D-instantons 

x 'ij — x 'i x j ^3.5^ 
with the momentum carried by the propagator in a corresponding QFT. 
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x k 




v i3 



Xjk Xj Xfc 



(a) 



(b) 



Figure 6: (a) D-instanton configuration and (b) corresponding momentum flow 

This identification would be incocnsistent if there was an interaction vertex which breaks 
their conservation law. In fact, any interaction vertex shown in Fig. || has following properties: 
(i) edges of the propagators are glued together in a definite cyclic order around the vertex so 
that successive edges share a common D-instanton label, (ii) all incoming momenta are given as 
the differences of XiS associated to successive boundaries. These two facts guarantee the sum 
of the incoming momenta is automatically zero^j For example, consider an interaction vertex 
among three D-instantons depicted in Fig. || (a). Since the incoming momenta are defined as 
Fig. |6] (b), their sum X{j + Xjk + Xki vanishes. It is now straightforward to check that Feynman 
rules given in Figs, [l] and [2] exactly coincide with those of a usual Yang-Mills theory. 

3.3 From sums to integrals 

In an ordinary perturbation method, we need to integrate over interaction positions in d- 
dimensional spacetime, which results in loop integrals. In the matrix perturbation theory, we 
need to sum over intermediate D-instanton positions. The correspondence between these two 
implies the equivalence of perturbation theory, because we have just seen that each fatgraph has 
the same factors both in Matrix and QFT pictures. 

As discussed in section 3.1, D-instanton configuration along flat M. d can reproduce d-dimensional 
momentum space. More concretely, let us assume D-instantons fill uniformly a d-dimensional 
hyperplane 



H := {(x 1 ,... , x d , 0, . . . ,0) £ M D } 



(3.6) 



3 The momentum conservation is a direct consequence of U(1) N gauge invariance. It is, however, nontrivial 
whether one can represent the momenta as differences successively. 
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with density p (i.e. there are p D-instantons per unit d-dimensional volume). Then, we can 
replace the sum over brane positions Ylk by the loop integral p f d d p. 

For example, in the case of the two point function {X^-X v -^ of section we can choose the 
point Xi as the origin of the momentum space. Then the translation dictionary reads 

Matrix QFT 

==>■ q (external momentum) 

Xjk P (loop momentum) (3 y\ 

Xik => Q+P 

: P [d d p... 

k 

Then eqs. (|2.11|) , (|1|), ( ^14]) can be written as 



$1+ 



1 



{q + p) 2 ■ p 2 

X {5^(q - p) X + 5 KX (2q + pf + 5 X ^{-q - 2p) K } 



x{S KU (p - q) A + 5 uX (p + 2q) K + 5 XK (-2p - q) v } 



( b ) + ( C ) = --IT I ^ + ) C 1 - d ) V 



(d) + (e) = -¥/^f^ 
9 7 {q+p) z p z 

(f)+(g) = -^J dp ^T>y? (3 ' 8) 

which look more familiar as those in standard QFT textbooks. 
3.4 Non-Abelian gauge symmetry in Matrix model 

Now we come back to the problem how we can incorporate the non-Abelian gauge symmetry 
of ci-dimensional SYM theory starting from Matrix models. Actually, without this non-Abelian 
structure, we cannot explain why the cyclic order of momenta is important in the amplitudes 

To achieve this, we need to consider the coincident D-branes as in Suppose we want 
to realize U(n) gauge symmetry. We need to put n D-instantons at the same point in M. D . 
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Hereafter the word "cluster" will be used to designate the n coincident D-instantons. The N 
D-instantons are thus grouped into M = N/n clusters. 

We choose a background in which cluster r is located at x r = (xj., ... , xf, 0, . . . ,0) £ R D : 



x£l 







































\ 



(/i = 1, . . . ,d), 



(3.9) 



V 

X m =0 (m = d + l,... ,D). 

Here 1 and denote unit and zero matrix of size n, respectively. 

In this background, U(N) gauge symmetry is broken down to U(n) M generated by 



U 



V 







































Umm 



\ 



(3.10) 



/ 



As for fluctuations, it is useful to divide N x N matrix into the blocks of size nxnas 



X 11 



X 11 
^11 






A 1M 


















vfJ, 






MM 



\ 



= 1, ... ,d), 



X 7 ' 



vm 
A ll 






ym 
A 1M 






ym 












vm 
A M1 






vm 



(m = d+l,... ,D), 



(3.11) 



•011 












As 












IpMl 









I 



4 We use r,s, . . . 6 {1, . . . , M} for cluster indices. i,j,... £ {1, . . . , N} are reserved for D-instanton labels. 
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Figure 7: D-instanton configuration corresponding to d-dimensional SYM. Each cluster is de- 
noted by a blob. The clusters are distributed along M. d with average spacing b. 



Each n x n block Xr S , VVs transforms as a bi- fundamental representation of U(n) r x U(n) s 
subgroup of U(n) M . As we will see shortly (cf. eq. ( |3.14| ) below), the fluctuations Xr S (fx = 
1, . . . ,d) tangent to H will be identified with the U(n) gauge field whereas those in the transverse 
direction X™ (m = d + 1, . . . , D) will play the role of Higgs fields. Similarly, D-dimensional 
spinor ip rs will be identified with their super partners. 

In this notation, U(n) M invariant correlation functions are something like 

(trX% r2 X% rs ...X% n ). (3.12) 



The correspondence given in sections |3.1|, 3.2 and 3.3 is true to this non-Abelian case, provided 



"D-instantons" are now replaced by "clusters." 

3.5 (i-dimensional super Yang-Mills action from Matrix model 

So far, we have studied how a gauge theory is embedded into Matrix theory, through the corre- 
spondence of correlation functions. 

To complete our analysis and to extract further intuition, it will be useful to rewrite the 
original Matrix model action into that of QFT: d-dimensional reduction of SYM in D dimensions. 
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In the notation given in ( |3.9| ) and ( |3.11|) , the action ( |2 . 1| ) reads 

S = ——y-'y^tT^-(x% s Xrs ~ x rs-^rs + [X^, X V ] rs ) (x^ r Xg r — X U gr X^ r + \X^ ,X V \ 



r.s 



+\{x% s X™ + [X»,X™) rs )(x» r X™ + [X", 

+i[X m ,X n ] rs [X m ,X n ], r (3.13) 

+^ sr r^x^ rs + [x^,p] rs ) 

+^ sr Y m [X m ^\ rs ) 

where the trace "tr" is taken over n x n matrix indices. 

Let us approximate the sum by the d-dimensional integral, which looks more like a QFT. As 



in section 3.4, we assume the clusters {x r } are uniformly distributed on d-dimensional hyperplane 
H with a constant cluster density p' . Since each cluster consists of n D-instantons, p' is related 
to the D-instanton density p as np' = p. 

Renaming the n x n matrix valued fields as 

X& =► X»{p) = A»{p) 

Xfs X m (p) = <S> m (p) (3.14) 

and the continuum approximation similar to (|3.7|), 



(3.15) 



the action fl3,13| ) reads 
S = -^± p i jd d p 



tr{i(^X» -p v X»(p) + [X»,X»](p))(-p»X»(-p) +p v X»(-p) + [X^X"]{-p)) 
+ \{p»X m {p) + [X»,X m ](p))(-p»X m (-p) + [X^X m ]{-p)) 
+ - A [X m ,X n ]{p)[X m ,X n }{-p) 
+ ^(-p)T^p^(p) + [X^ij}(p)) 



+ 2 



l -^-p)T m [X m ^]{p)). 

(3.16) 
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This is a momentum space representation of the super Yang-Mills action in (^-dimensions. 
Using formula like 

-i [ d d x(d^A u (x))e- ip - x 



x»{p)X v {p) 
[X^X v ](p) 



p'{2-K) d / d d x[A»(x),A u (x)]e- ip - x , 



(52 r jj Jfo^\d 
2 P 

g 2 

x tr 



d d x 



inverse Fourier transform of ( p. 16|) gives 
5 

i {(d^A^x) - d u A»{x)) + ip'(2ir) d [A^(x),A u (x)}y 
+ ±{(d^ m (x)+ip f (2Tr) d [A»(x),$ m (x)]y 
-i{p'(2vr) d [$ m (x),$ n (x)]} 2 

' $(x)T»{d^{x) +ip'(27r) d [A >1 (x),iP(x)}} 



Rescaling the fields as 



- 1 -{i,{x)T m p'(2^ m {x)Mx)]} 



p'(2Tr) d A^(x) -» A»(x) 
p'(2ir) d $ m (x) -» $ m (x) 

and defining the (i-dimensional coupling constant by 



(3.17) 



(3.18) 



(3.19) 



(3.20) 



we are finally lead to a familiar form of (i-dimensional SYM coupled with adjoint matters :Q 

S d = J d d x\tl^{F^) 2 + i(ZV$m)2 _ ^[$™ $ n]2 
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+^(zyv)-^rn* m ,V>] 



(3.21) 



5 We neglected the factor (yj 1). This reflects the fact that overall shift of D-instanton configuration results 
in the same momentum configuration in SYM picture. To make the mapping one-to-one, we need to specify the 



origin of momentum space as discussed toward the end of section 3.f . (Admittedly, the replacement (3.7) is some- 
what misleading.) Other possibility is to introduce an additional gauge symmetry to constrain the off-diagonal 
components as X r<s = X r+ k, s +k (Vfc £ Z). This gauge symmetry kills the ambiguity of overall shift, but the 
clusters need to be arranged periodically on a lattice. The latter method is equivalent to the S 1 compactification 
proposed by W. Taylor BL4J. 
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Here, the standard covariant derivative for the adjoint matter D^$ m = d^Q m + i[A^,^ m ] and 
the field strength F^ v = d»A u - d v A» + i[A^, A v \ is used. 
Our procedure can be schematically summarized as 




(3.22) 



Now we make a few comments on the hermiticity of the fields. If we naively interpreted the 
brane configuration as the discretization in the spatial coordinates (a la lattice gauge theory) 
rather than the momentum coordinates, then annxn block Xrs must be interpreted as a U (n) 
connection (or parallel transport) matrix connecting two points x r and x s . Then it must be anti- 
hermitian (or unitary) matrix. But, since Xf. s is just a part of a much larger N x N hermitian 
matrix, there is no a priori reason why it should take such special forms. It might be possible to 
devise a mechanism to put such a constraint, it would lead to additional complication. Actually, 
the hermiticity as a Ax N matrix leads to (Xr S y = Xs r . This is consistent with our identification 
Xrs with A fl (p), because the hermiticity of the gauge field A*(x) = j (^) d (A^(p)e ip - x + /i.e.) 
leads to the condition (A fl (p))^ = A^(—p). 

Note also that the diagonal parts X£ of N x N matrix field correspond to the zero-momentum, 
Cartan component of U(N) gauge field. 



4 Divergences in large N matrix theory 

Although the Matrix model is defined for arbitrary value of N, interesting physics is believed to 
emerge from the large A limit. But there seems to be no general argument or rule concerning 
what kind of large N limit should be taken. 



It is well known that in the so-called 't-Hooft limit 12] N — > oo with gy^N being fixed, 
U{N) Yang- Mills theory simplifies drastically, and some exact analysis, say Borel summability 
became possible. Eguchi and Kawai fi~5[| have argued that four dimensional large N gauge theory 
can be replaced by Matrix models in zero dimensions. In a sense we are studying the reverse 
process of the Eguchi-Kawai reduction from a different viewpoint. 

In type IIB matrix models, g appears as the overall factor in the action and thus there is 
no nontrivial critical point for g. Yet, it is still controversial fl~6fl how large N limit should be 
taken; g fixed? g 2 N fixed? or what else? 
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Figure 8: D-instanton configuration corresponding to (i-dimensional SYM. N D-instantons are 
distributed in a (i-dimensional ball with radius A within M. d 

In the case of SYM theory embedded in the matrix perturbation theory, one is faced with 
an additional complication. Increasing N allows two different interpretation: (i) larger gauge 
symmetry (ii) larger cutoff (i.e. more degrees of freedom). In the usual QFT case, they are 
clearly separated and never mixed. Here, however, there is a crosstalk between the two. 

In the perturbative computation of (i-dimensional quantum field theory, divergences arise 
from the integrals over loop momenta f d d p. We can regularize the integrals by introducing a 
UV cutoff A. Suppose the amplitude associated with a Feynman diagram T diverges as 
in the continuum limit A — > oo. Basically, D (T) can be determined from the superficial degree 
of divergence of T or its subgraphs. Renormalizability is the property that all divergence can 
be removed if A — * oo limit is taken not keeping g fixed but adjusting g so as to maintain a 
certain functional relation H(g,A/[i) = 0. It is of course a very nontrivial problem to find the 
explicit form of renormalization trajectory H(g,A/(i) = 0, but at least perturbatively, it can be 
determined by carefully analyzing the structure of divergences in Feynman diagrams. 

On the other hand, in the matrix perturbation theory we are working with, only source of 

divergence is sending N to infinity.^ In the spirit of correspondence between Matrix model and 
6 Of course, there are combinatorial divergences due to the infinitely many Feynman graphs. But this is common 
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5 l/2 




core size 


the minimum distance between D-instantons 


a 




D-instanton spacing 


average distance to the nearest D-instantons 






renormalization scale 


typical momentum scale of the external lines 


A 




cutoff 


distance to the farthest D-brane 



Table 1: Various scales in D-instanton distribution 

QFT, continuum limit A — > oo should be related to the large N limit. In other words, adding 
more and more D-instantons on the outskirts of the D-instanton cluster should be equivalent 
to increasing A. In this picture, the UV cutoff A is nothing but the distance to the farthest 
D-instanton (see Fig. [?].), which is natural from UV/IR correspondence 12] or spacetime uncer- 



tainty JTT|. 

The analysis of large N behavior can be complicated because there is no unique way to add 
extra D-instantons; the relation between the two limits N — > oo and A — > oo is highly dependent 
on the strategy of putting new D-instantons. 

In this section, we will take a phenomenological approach to clarify the relation between 
large N limit and continuum limit A — > oo. N D-instantons are assumed to be concentrated 



along d-dimensional hyperplane with uniform density p as in section 3J3. But here, since N is 
finite, the radius A of D-instanton cluster is also finite. 



4.1 Scaling laws 

So far we have encountered various length scales. The core size g 1 ^ 2 , the D-instanton spacing 

a, the scale of external momenta fi. For finite N, the distance A to the farthest D-instanton 

will also play an important role. We will study the scaling laws for these length scales. For the 

readers' convenience they are listed in Table 1. Although the word "length" will be frequently 

used, they represent momentum scales in the QFT picture as we argued in the previous section. 

We are interested in how they should be varied as N tends to infinity. In this paper, we 

assume a simple power law scaling and try to draw some bounds on the exponents from physically 

reasonable assumptions. 

Since the quantities in Table [l] are all dimensionful while N is dimensionless, one must decide 
to Matrix and to both theories and is not discussed in this paper. 
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b 


e 


-0 


cluster spacing 


P 




d.0 


D-instanton density 


pf 


~N e 


d0 
A 1 " 9 


cluster density 


9l 


~ N -" 


d{9-w) 

r-u A 


YM coupling in (f-dimensional QFT 



Table 2: Other quantities with nontrivial scaling laws 



which is kept fixed in the large N limit. For this purpose we choose //, the momentum scale 
carried by the external lines in the QFT picture. In other words, all "lengths" discussed in this 
section are measured in the unit of p,. For example, we set p ~ Xij when we compute (X^X^). 
The other three quantities g 1//2 ,a and A are assumed to scale with some power of which is 
specified by three independent exponents d, 9 and u as in Table |l|. 

In addition to these basic "length" scales, there are some other quantities of interest, with 
nontrivial N dependence: 

• D-instanton density p 

For finite N, D-instantons are assumed to be distributed uniformly with density p within 
a d-dimensional ball of radius A. Thus we have 

pA d ~ N (4.1) 

Since A ~ N < - 1 ^/ d , (|4~l| ) fixes the scaling of the p as 

p~N 9 . (4.2) 

• cluster spacing b 

The D-instanton spacing a and cluster size b are related via na d = b d . Since we fix the 
rank n of the gauge group, a and b will have the same scaling behavior, b ~ a ~ N~ e / d . 

• cluster density p' 

By the same token, the cluster density p' , related to D-instanton density p via np' = p will 
have the same scaling as p. 

p' ~ N 9 (4.3) 
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• YM coupling in d-dimensions gd 

YM coupling in d-dimensions is related to g and p' via ( |3.20| ). Thus it will have a 
nontrivial ./V dependence: 

g 2 ~ p'g 2 ~ iV e - (4.4) 

We summarize the result in Table |2[ 

4.2 Physical bounds on scaling exponents 

We have seen that a d-dimensional SYM theory emerges from the off-diagonal dynamics of the 
large N Matrix model. In order to prove the claim, we need at least to show such a large 
N limit is indeed possible — precisely specifying how to arrange the background D-instanton 
configuration as iV tends to infinity. It may be difficult to do this rigorously. We will content 
ourselves with obtaining some inequalities among the exponents 9, u, d introduced in section 



44], so that there occurs no apparent inconsistency in QFT side. This would help us applying 
Matrix theory to more realistic situations in the future. 

We will consider several physically reasonable assumptions, but we do not intend to claim 
that following conditions are all necessary or sufficient. 

(i) Interpretation as d-dimensional QFT required replacing the discrete sum by (i-dimensional 
integral. This coarse graining can be justified only when a ~ b <C p. This is true in the large N 
limit if 

6 > 0. (4.5) 

(ii) From the d-dimensional point of view, cutoff scale A must tend to infinity. Thus we 
have 

e < l. (4.6) 

(iii) As we saw in section ^, the perturbation theory is essentially the expansion in g/a 2 . 
Thus it is valid if a 3> g 1 ^ 2 is satisfied. This remains to be true in large limit if 

dcu>49. (4.7) 
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(iv) Actually, the same bound can be obtained from a different viewpoint. From exact 



results for matrix integrals [17, 18, 19], it is reasonable to assume there is a pairwise repulsive 
potential among D-instantons due to entropy factor. This could be effectively treated [|| as 
each D-instanton has a core size of order g 1 / 2 . This implies a > g 1 ^ 2 ■ Sending N to infinity, we 
obtain the inequality (fO|). 

(v) To construct U (n) gauge theory, n D-instantons are put on the same point (see section 
3.5). But this assumption might be too strong; it is possible that the n D-instantons can 



disperse in the cluster of size b but are still grouped via (slightly broken) U(n) gauge action 
from QFT point of view.Q Of course, the dispersion size b must be sufficiently smaller than the 
renormalization scale, 

b < ii. (4.8) 

This corresponds to the minimum momentum resolution seen by U(n) Yang-Mills theory. The 
condition (4.8) leads to the bound, 



(vi) As far as a tree level amplitude or the form of Yang- Mills action is concerned, the argu- 
ment given in (v) is sufficient. But if quantum effects are taken into account, it is another story. 
Just like the anomaly from one loop, large N divergence from loop integrals may overwhelm the 
b/n ~ N~° suppression discussed in (v) and may yield non-negligible effects. 

Let us estimate the effect of dispersing D-instantons using the one-loop two point function as 
an example. As a function of D-instanton configuration {xk}, the most divergent contribution 
is roughly given by 



,l-loo P rr_ ,i 9 2 1 9 2 P f A d d p 



•A p [{ a; fc}] ■ ■■> / ■-> > i ■-> > 
V- t x tk V J P + V 



~ ll^- 2 



(4.9) 



7 This claim is not so strange as it sounds. In Nature, non-Abelian symmetry is exact in UV regime but hidden 
in IR regime through confinement or Higgs mechanism. In our context, the D-instantons within a cluster look 
almost coincident in much larger scale y,. 
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If the D-instanton positions {xk} have a dispersion of order b, A op [{xk}] will change as 



The last expression scales as N to the (9 — co — % + (d — 3)— ^-)-th power. Thus, it is negligible 
only if 



This is the condition when the operators associated with the slight shift of D-instantons are 
irrelevant in the sense of large N renormalization group. 

(vii) Let us study the problem of renormalizability of d-dimensional theory i.e. whether or 
not as N — > oo only a finite number of amplitudes superficially diverge. 

From the Wilsonian point of view, renormalizability is not a necessary condition for QFT, 
but a consequence of the renormalization procedure. But it is of some interest in presenting the 
analysis since in the context of Matrix model, situation is rather complicated 

For a given Feynman graph, the superficial degree of divergence is usually determined from 
the number of loops and propagators. The coupling constants just count the number of vertices 
and stay fixed when the cutoff is sent to infinity. As we all know, d = 4 is the critical dimension 
for gauge theories. 

But here, we are talking about the divergence when N tends to infinity. Recall not only A 
but also gd changes as a function of N, i.e. g\ ~ N d ~ u . Thus usual definition of the superficial 
degree of divergence does not work. 

In a sense, we are studying a generalized large N limit in which g^N u ~ e is kept fixed. Thus 
the standard QFT results should follow if we restrict to 9 = u, whereas 't Hooft limit would 
correspond to another special case, 6 + 1 = u. 

What is the new rule for the superficial degree of divergence? Note that g 2 always come in 
pair with a sum over D-instantons. From the substitution 




(4.10) 



d(l -u) < 3- 29 



(4.11) 




k 



(4.12) 



8 



For example, it is not clear "What is IR limit?' 
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extra (i-dimensional loop integral is always associated with the factor g\. Each loop contributes 
A d while the coupling gives gj ~ N 9 ^ ~ \d(e-u,)/(i-e) It 

is easy to convince oneself that the 
net effect is to replace the spacetime dimension d by an effective dimension 

d(8 — u) ,1 — u , , , 

d es^ d + \^r = d T^e- (413) 

Thus we have a new criteria about large iV renormalizability as follows 

d eff < 4 ^ d(l - u) < 4(1-0). (4.14) 

As promised, 6 = uj recovers the standard result. Note that if u = 1, d e jf = for any d, 0. This 
corresponds to the well known fact that planar limit of the + O-dimensional Matrix model 
absolutely converge. 

(viii) As for the error in replacing sums by integrals, analysis in (vii) can be generalized to 
an arbitrary Feynman graph T. Suppose we know the amplitude diverges as 

A(T) ~i\P< r \ 

including N dependence of g\. Then, the approximation can be justified if 

5A ~ 4iV 7(r) ~ JSP^-a « i. 
A 

Thus the error is negligible for graphs with sufficiently low degree of divergence: 

7 (r) < - d . (4.15) 
It may be useful to introduce effective superficial degree of divergence, V e ff(T) defined through^] 



or equivalently 



Then, ( gig) can be expressed as 



V^T) < ^ 6 > 1 - ^ (4.16) 



T> e ff(F) coincides with usual superficial degree of divergence T>(T) if gd is independent of N. 

26 



Suppose (i-dimensional theory is renormalizable in the sense that D e ff(T) has a T-independent 
upper bound T> max . Then by choosing 

B>1 1 



V 



the error can be neglected for all Feynman integrals. In particular, d = 4 SYM case (T> max = 2) 
leaves finite window 1 > 9 > tt. 

If ci-dimensional theory is non-renormalizable, T> e jj grows with the number L of loops. Thus 
from ( |4.16| ), for a given 9, there is a maximal number of loops L max beyond which the approxi- 
mation fails. But L max tends to infinity if we approach 't Hooft limit, 9 — » 1. 

4.3 Possible interpretation of exponents 

So far we have chosen a particular D-instanton configuration depicted in Fig || the dimension 
d in which QFT lives is determined by the number of flat directions. But this is clearly a very 
special configuration from D-dimensional viewpoint. Can we relax the assumption? 

In Wilson's approach to renormalization group, one can study the origin of ultraviolet di- 
vergences by isolating the dependence of the functional integral on the short distance degrees of 
freedom of the field. In Matrix approach, short distance degrees of freedom correspond to the 
long distance D-instantons. The number SN of D-instantons contained in the momentum shell 
A < \p\ < A + 5A is given by 

6N oc A d ^5A as A -» oo. 

In fact, as far as the loop divergence is concerned, dimension d will appear only through this 
relation. 

Consider for example a more generic D-instanton configuration, Fig |9[ Let N(A) be the 
number of D-instantons within a L>-dimensional ball of radius A. Then the space-time dimension 
can be "defined" as the rate of growth of N: 

rf = ^£^>. (4,7) 
olog A 



For a uniform configuration like Fig ^, definition ( [4.17 ) gives the number of flat directions. Note 



that the new definition ( 4.17 ) make sense for non-integer dimension d, which may be useful to 



visualize the meaning of dimensional regularization. 
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A + 5A 




Figure 9: General D-instanton configuration. The number 5N of D-instantons contained in the 
shell A < \p\ < A + 5 A is proportional to A d_1 . 

In order to enumerate the physical degrees of freedom in a field theory, one needs to put the 
system into a finite box of volume V. The number of states is given by the available phase space 
volume. 

In order to realize a continuum field theory in infinite spacetime, one needs to take two limits: 

V — > oo (Large volume limit), 

A — ► oo (Continuum limit). 
Singularities associated with the former and latter are usually called IR and UV divergences, 
respectively. In standard textbooks on QFT, V — > oo limit is taken first so that Feynman rules 
simplify in the momentum space. Subsequently, A — ► oo limit is carefully investigated. This 
asymmetry between the two limits is due to the well known fact that the translational invariance 
in momentum space is actually broken by the hierarchical structure. 

In the Matrix model, however, all limiting processes are "unified" into a single large N limit. 



Comparing ( 3.22 ) and ( 4.18 ), we can say that we have investigated in section O all possible 
limits 

V ~ (J ~ N 9 , A ~ N^-°^ d (4.19) 

to get a continuum field theory. 
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5 Discussions 



In this paper, we have studied the dynamics of large N Matrix models through the quantum 
fluctuations in a fixed D-instanton background. In particular, we have explicitly shown the 
correspondence of perturbation theories between the usual QFT and Matrix perturbation theory. 
The correspondence is exact if relative D-instanton positions are interpreted as momenta in QFT 
picture. 

One might think that this is a kind of triviality. Indeed, Matrix model action is originated 
from the Yang-Mills action by dimensional reduction. It is no wonder Yang-Mills theory can 
be recovered from the Matrix model. However, since dimensional reduction is simply throwing 
away spacetime coordinate dependence, the reverse procedure would be just re-introducing x 
dependence to the matrix fields. But contrary to this naive expectation, the momentum space 
picture emerges first and coordinate picture is recovered only after Fourier transformation. 

The correspondence exploited in this paper can be regarded as a "dual" version of Eguchi- 
Kawai reduction [pi]]. The original suggestion by Eguchi and Kawai is valid only at strong 
coupling pQ l, whereas we have shown the equivalence in a weak coupling regime. 

Just like lattice gauge theories, Matrix model provides us with a natural gauge invariant reg- 
ularization. But "Matrix regularization" has two important features. First, the quantum fields 
are discretized in the momentum space rather than ordinary space, and the hierarchical structure 
inherent to QFT can be understood in a geometrical fashion. Second, matrix regularization can, 
in principle, be "generally covariant" if the sum over all background configurations is taken into 
account. A permutation of D-instantons is a discrete analogue of coordinate reparametrization. 

Matrix models pack too much degrees of freedom into a few matrices. As is often the case, 
this obscures the meaning of large N limit. Furthermore, in a theory with T-duality it is difficult 
to make distinction between IR and UV limits. The limits explicitly depend on the effective 
dynamics we are talking about. At any rate, it is obvious in Matrix theory that universal 
behavior is expected only in large N limit. 

We initiated a preliminary study of what class of large N limit is possible in order to repro- 
duce a QFT. Key idea is to classify the degree of divergence in terms of N, the only source of 
divergence in Matrix theory. It is now possible to interpret renormalization group a la Brezin 



and Zinn- Justin f2U 22 , \23 in terms of the usual renormalization of Yang-Mills theory. We hope 
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to report on this elsewhere. 

One might well be puzzled by the interpretation of the Matrix dynamics as d-dimensional 
Yang- Mills theory If the background D-instantons are distributed really uniformly the rank of 
the gauge group would be just a matter of choice because it depends how we cut the D-instanton 
gas into pieces. Let p be the D-instanton density. Suppose we decide to call D-instantons inside 
a ti-dimensional hypercube of size b as a cluster. Then we have 

n = pb d , (rank of the gauge group) 

p' = p/n = b~ d , (cluster density) 

g j = (2ir) d p'g 2 = (2ir) d g 2 b- d . (Yang-Mills coupling) 

Note that b dependence cancels in the 't Hooft coupling A = g\n = (2n) d g 2 p. Therefore, matrix 
perturbation theory suggests that any universal property of U (n) gauge theory with adjoint 
matters should depend, not separately on g^ or n, but on 't Hooft coupling A = g\n, at least for 
sufficiently small A. 

At present, we do not know whether this is generally true or not, but the following evidence 
should be taken seriously. Consider a renormalization group beta function for d = 4 Yang-Mills 
coupling, 

P(94) = -A)<?4 - 0104 • 

In a U(n) gauge theory with adjoint matters ((^(G) = n), the coefficients are given by (3q = 
cqu, 0i = c\n 2 , . . . with cq, c\, . . . depending only on the matter contents. Thus we have 

(3(g A ) = -c ngl - a n 2 g\ . 

This can be rewritten as 

/3(A) = -2(c A 3 + ci A 5 + •••)• 
The right-hand side is a function of A only, in accordance with our expectation. 

This work is supported in part by the Grant-in- Aid for Scientific Research on Priority Area 
707 "Supersymmetry and Unified Theory of Elementary Particles", Japan Ministry of Education. 
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